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Applying random matrix theory to quantum transport in chaotic cavities, we develop a novel approach to 
computation of the moments of the conductance and shot-noise (including their joint moments) of arbitrary 
order and at any number of open channels. The method is based on the Selberg integral theory combined with 
the theory of symmetric functions and is applicable equally well for systems with and without time-reversal 
symmetry. We also compute higher-order cumulants and perform their detailed analysis. In particular, we 
establish an explicit form of the leading asymptotic of the cumulants in the limit of the large channel numbers. 
We derive further a general Pfaffian representation for the corresponding distribution functions. The Edgeworth 
expansion based on the first four cumulants is found to reproduce fairly accurately the distribution functions in 
the bulk even for a small number of channels. As the latter increases, the distributions become Gaussian-like in 
the bulk but are always characterized by a power-law dependence near their edges of support. Such asymptotics 
are determined exactly up to linear order in distances from the edges, including the corresponding constants. 
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I. INTRODUCTION 

Quantum transport in mesoscopic systems has been a sub- 
ject of an intense study during the last decade, see reviews'-^. 
Traditionally, the focus of interest has been on statistical prop- 
erties of the conductance, g, and shot-noise, p. For noninter- 
acting electrons, the Landauer-Biittiker scattering formalism^ 
relates these quantities (expressed in their natural units) to the 
so-called transmission eigenvalues Ti of the conductor: 

n n 

ff^^T,, p = Y,T^{l-T,). (1) 

i=l i=l 

Here, n = min(A^i, iV2), where iVi.2 are the number of prop- 
agating modes (channels) in the two attached leads. The T, are 
the eigenvalues of the matrix tt^, with t being A^i x N2 matrix 
of transmission amplitudes from the entrance to exit channels. 
They are mutually correlated random numbers, < < 1, 
whose distribution depends on the type of the conductor. 

Below, we consider chaotic cavities (open quantum dots). 
In this case, random matrix theory (RMT) has proved to 
be successful in describing universal fluctuations in transport 
through such systems.'"* Within this RMT approach, the joint 
probability density function (JPDF), Vp{T), of the transmis- 
sion eigenvalues is induced by the random scattering matrix 
drawn from one of Dyson's circular ensembles^, according 
to the global symmetries present in the system. The exact 
expression for this JPDF is knowni^ to have the following 
simple form: 

n 

Vp{T)^M^^\HT)f\{Tr\ (2) 

1=1 

where a = |(|A^i - N^] + 1) and A(r) = Y{,^j{T, - T,) 
denotes the Vandermonde determinant. The Dyson's sym- 
metry index (3 depends on the presence {j3 ~ 1) or absence 



(/3 ~ 2) of time-reversal symmetry or that of spin-flip sym- 
metry {(3 — 4) in the system, thus distinguishing between 
the three canonical RMT ensembles (orthogonal, unitary or 
symplectic, respectively).^ The normalization constant Mp is 
given by 

_'fr r(i + |(i + j))r(a + |j) r(i + |j) 

/i r(l + f)r(l + a + f(n + j-l)) 

and assures that expression (|2| is a probability density. It is 
known generally for discrete positive n and continuous a and 
13 as the Selberg integral. 

Presently, there is a substantial progress in the understand- 
ing of statistics of various transport observables in chaotic 
cavities that is due to the recent developments of new ana- 
lytical methods in the theory. Among them, the Selberg in- 
tegral theory plays a special and important role, see the re- 
cent review^ for its current status. In the present context, it 
has been initially applied by two of us2 to find the average 
value of shot-noise and hence the Fano factor exactly. This ap- 
proach has then been developed further to study full counting 
statistics of charge transfer'" as well as to obtain exact explicit 
expressions for the shot-noise variance and for the skewness 
and kurtosis of the charge and conductance distributions"^!^ 
(see also Refs.'3i'^-'^-'^i'A"^i'V" for other RMT results on the 
relevant statistics in chaotic cavities). Since j3 enters the Sel- 
berg integral as a continuous parameter, this method allows 
us to treat all the three ensembles on equal footing, thus giv- 
ing a powerful alternative to diagrammatic'^''^ or orthogo- 
nal polynomialii2£ approaches, especially when the channel 
numbers are small. 

A completely alternative treatment has been recently under- 
taken within the semiclassical approach^ii^^ which represents 
quantum transport in terms of classical trajectories connect- 
ing the leads. By constructing asymptotic semiclassical ex- 
pansions for transport observables, this approach successfully 
accounts for both system-specific and universal (RMT) fea- 
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tures, see Ref. for a review. 

The case of (3=2 (broken time-reversal symmetry) is known 
for several reasons to be the special one in RMT. For the prob- 
lem in question, further progress in this case has been made 
very recently along the following two directions. Novaes2i 
combined the Selberg integral with facts from the theory of 
symmetric functions to compute non-perturbatively moments 
of both the transmitted charge and conductance but not those 
of shot-noise. Alternatively, Osipov and Kanzieper^^ com- 
bined the theory of integrable systems with RMT, as given 
by dU, bringing out an effective formalism for calculating the 
cumulants of the conductance and shot-noise. However, the 
relevant consideration for the systems with preserved time- 
reversal symmetry, /3=1, is still lacking. 

The distribution functions of the conductance and shot- 
noise are also studied quite intensively on their own. How- 
ever, no explicit expressions have been reported so far ex- 
cept for a few cases, namely, for the conductance distribu- 
tion at Ni^2 = 1, 2^'^''^'^^ and for the shot-noise distribution 
at = Asymptotic analysis of the both distributions 
at '3> 1 has been performed very recently in Ref;^, see 
also Ref.^^. To the best of our knowledge, no general results 
valid at arbitrary A^i 2 and /3 7^ 2 are available thus far Mean- 
while, the conductance distribution with dephasing^" has been 
directly measured in open quantum dots^' and in microwave 
billiards. The shot-noise power in chaotic cavities has been 
recently studied experimentally.^^ Counting electrons in quan- 
tum dots is also experimentally accessible..^ All this provides 
an additional motivation for the present study. 

In this work we explore further the direction along the lines 
of Novaes's work^'* and develop a systematic approach for 
computing the moments of linear statistics in transmission 
eigenvalues for the systems with both preserved and broken 
time-reversal symmetry. This approach yields the moments 
of the conductance and shot-noise of arbitrary order, includ- 
ing their joint moments and cumulants. In the next section 
we present the detailed exposition of the method used, includ- 
ing the relevant facts from the theory of symmetric functions. 
This method is then applied in Sec. |III]to derive expressions 
for the moments and cumulants of the conductance and shot- 
noise in a closed form. Sections |IV] and [V] complement this 
study by investigating the corresponding distribution func- 
tions and their asymptotic behaviour Our main findings are 
summarized and discussed in the concluding section [Vll 



II. THE METHOD 

The method is based on expanding powers of the conduc- 
tance or shot-noise (or any other linear statistic) in Schur func- 
tions s\{T). These functions are symmetric polynomials in 
the transmission eigenvalues T = {Ti, . . . , T^} indexed by 
partitions A. In the group representation theory the Schur 
functions are the irreducible characters of the unitary group 
and hence are orthogonal. This orthogonality is quite useful 
since it means that the coefficients in Schur function expan- 
sions are just "Fourier coefficients", and, hence, can be found 
by integration over the unitary group. It gives an efficient 



way of calculating the expansion coefficients explicitly, the 
fact that we exploit in our approach. The Schur functions can 
be then averaged over the JPDF (|2]i, 

/n 
d[T]sx{T)VpiT) , d[T] = n dT, , (4) 
i=i 

with the help of integration formulas due to Hua-^^. The Schur 
function expansions and Hua's integration formulas provide 
us with the necessary ingredients to compute all moments (or 
cumulants) of the conductance and/or shot-noise, see Sec. |III] 
for the detailed analysis. 

In this section, we first give a brief summary of the required 
facts about partitions and Schur functions, then develop the 
systematic way of performing the expansion over Schur func- 
tions and finally determine Schur function averages. 



A. Partitions and Schur functions 

A partition is a finite sequence A = (Ai, A2, . . . , Am) of 
non-negative integers (called parts) in decreasing order Ai > 
A2 > ■ ■ • > Am > 0. The weight of a partition, |A|, is the 
sum of its parts, |A| = J^j ^j^ ^"d the length, 1{X), is the 
number of its non-zero parts. No distinction is made between 
partitions which differ only by the number of zero parts. Dif- 
ferent partitions of weight m represent different ways to write 
m as the sum of positive integers and can be graphically vi- 
sualized through the Young diagrams. For example, one has 
only one partition A = (1) in the trivial case of m = 1; 
two partitions A = (2,0), (1, 1) for m = 2; three partitions 
A = (3, 0, 0), (2, 1, 0), (1, 1, 1) for m = 3, etc. 

For any partition A of length /(A) < n, one can define a 
symmetric polynomial sx in n variables xi, . . . ,Xn as fol- 
lows; 

det|:r^+"-^r 
SA(a;i, . . . ,a;„) = r^j (5) 

The denominator here is nothing else but the Vandermonde 
determinant A(x) = Y[i<ji^i ~ ^j)- divides the corre- 
sponding factor in the nominator, leaving the quotient as a ho- 
mogeneous polynomial in the Xj's of degree m = |A|. These 
polynomials s\ are called the Schur functions. For one-part 
partitions, A = (r), Schur functions are just the complete 
symmetric functions, S(r) (x) = hr, while for partitions which 
have no parts other than zero or one, A = (1, . . . , 1) = (l*"), 
the Schur functions S(ir) are the elementary symmetric func- 
tions e,- (x). This can be verified directly from (|5]). It should be 
noted that the Schur functions corresponding to the partitions 
of m form a basis in the space of homogeneous symmetric 
polynomials of degree m, so that any homogeneous symmet- 
ric polynomial can be written as a linear combination of Schur 
functions. 

The Schur functions of matrix argument that we shall use 
below are defined by the rhs in (|5]l evaluated at the eigenvalues 
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of the matrix. Taking as an example the nxn matrix T = tt^ 
of transmission probabiHties, one has 



where Ti , . . . , r„ are exactly the transmission eigenvalues 
that appear in (|2|. Although not apparent from this defini- 
tion, the Schur functions of matrix argument are polynomials 
in the matrix entries-'^ and, obviously, s\{T) = s\{XTX^^) 
for any non-degenerate matrix X. 



B. Schur function expansions 

In order to determine the moments of the conductance and 
shot noise along the lines explained above, one needs to ex- 
pand the powers of these quantities in Schur functions. To 
this end, it is more convenient to work with the corresponding 
generating functions e* ^ "^^ or e* ^ "^^ (i^^j ) . These functions 
belong to the general class of multiplicative symmetric func- 
tions, where the coefficients of the Schur function expansion 

3 A 



(6) 



can be determined explicitly provided that the function / is 
analytic in a neighborhood of |a;| — 1 in the complex x-plane. 



see, e.g.. Appendix in Ref. . Indeed, thinking of the 
of the eigenvalues of a unitary matrix U , one can write 



(7) 



The main advantage of going unitary is the orthogonaUty of 
Schur functions (d/i(J7) is the normalized Haar measure): 



U{n) 



(8) 



which is a fact from the theory of group representations. One 
now recognizes a "Fourier series" in (|7|l and, hence. 



.(/) _ 



U{n) 



dfi{U)F{U)sl{U). 



(9) 



The integral on the rhs in (|9|l is a standard one in RMT. 
To evaluate it, one first transforms it to the eigenvalues 
e'^^ , ■ • ■ , e'^" of the unitary matrix U. The corresponding Ja- 
cobian is |A(e*^)p, canceling the denominator in the Schur 
function s^(e'^). The resulting integral can then be evaluated 
with the help of the Andrejeff identity yielding 



= det 



In view of the analyticity one can abandon the restriction 
I a; I = 1. Writing f{x) = '^jTjX^, one brings the Schur 
function expansion (|6]l and ( fTOl ) to the following general form: 



n / +C30 



n( E xA = Y,cx{r)sx{x), (11a) 

i—1 ^j— — oo A 

cx{t) = dei{Tx^-k+i}l (lib) 



The summation here is over all partitions A of length n or less, 
including empty partition (0) for which sx = 1. Expansion 
(fTTl i was also obtained by Balantekin"*" by algebraic manipu- 
lations. 



C. Schur function averages 

Hua in his book^^ evaluated many useful matrix integrals. 
The following two are relevant in the context of our work: 

/ \ TT ^(■^' + l)r(Aj- + n-j + a) -pr . 



3 = 1 ^ ■' J J l<i<j<r 



and 



{Sx)f!=l = 



2""!ni<.<,-<„(A.-A,-i+j) 
Ill<^<,<ni^^ + + 2n + 2a - i ~ j) ' 



(12) 



(13) 



where the average {sx)p is over the JPDF Q, as in Eq. (|4|l. If 
A = (0) then sx = 1 and both integrals follow from the Sel- 
berg integral^. It should be noted that the rhs in ( fT2] i is exactly 
Selberg's expression and the rhs in (T3[ can be manipulated to 
the one obtained by Selberg with the help of the duplication 
formula for the Gamma function. For non-empty partitions A, 
the integral in ( fT2] i is a particular case of the Kadell-Kaneko- 
Yan generalization"*'-'*^ "^^^ of the Selberg integral. However, the 
integral in ( fT3] ) is different, as the Kadell-Kaneko-Yan gen- 
eralization of the Selberg integral for /3 = 1 involves zonal 
polynomials. 

Hua's identities are valid for arbitrary continuous a > 0. 
Specifying further to our case of a = §{\Ni - N2\ + 1) and 
N = Ni + A^2, we arrive after a simple algebra at 

= {N,~jy. {N,-jy. 

(Xj + N-jyr 

where we have introduced the coefficient^^ 



(14) 



cx 



m=i m+x,-jy 



(15) 



(10) and 



{S\)f3= 



Cx 



X 



n 



(A,-+jV2-j)! 

\ {N2-jy 



n iwri 



l<i<j<N2 



+ l + Xi+Xj-i-j 



(16) 



The symmetry between iVi and N2 is not apparent in ( fTSI ). 
One can rearrange the terms in the second product on the rhs 
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TABLE I: Schur function expansion {{"^ Tj)^)i3 = m! ^ ca {s\)i3 for (3 = 


1,2. 




m 


Partition A 


(sa)/3=2 


{s \) 3=1 /{S\) 13=2 


C\ 


1 


(1) 


N1N2 
N 


N 
N+1 


1 


2 


(2,0) 


JVi(JVi + l)JV2(]V2 + l) 
2N{N-\-l) 


N+1 
N+3 


1 

2 




(1,1) 


JVi(JVi-l)JV2(]V2-l) 

7\7 / AT ~\ \ 

Zl\ [I\ — i) 


N-1 

N+1 


1 
2 




(3,0,0) 


iVi(iVi+l){]Vi+2)JV2(]V2+l)(]V2+2) 
6JV(]V+l){]V+2) 


N+2 
N+5 


1 
6 


3 


(2,1,0) 


(iVi-l)JVi(JVi+l)(JV2-l)iV2(JV2+l) 
3{]V-1)JV(]V+1) 


N 
N+3 


1 
3 




(1,1,1) 


(JVi-2){]Vi-l)]Vi(JV2-2){]V2-l)JV2 
6{]V-2)(JV-1)]V 


N-2 

N+1 


1 
6 




(4,0,0,0) 


JVl(iVl + l){JVl+2){iVi+3)iV2(iV2+l){]V2+2)(JV2+3) 
24]V(A'+l)(JV+2)(Ar+3) 


N+3 
N+7 


1 

24 




(3,1,0,0) 


(JVi-l)JVi(]Vi+l){]Vi+2)(JV2-l)]V2(JV2+l)(iV2+2) 
8(A^ — l)iV(A^-|-l)(A^+2) 


N+1 

N+5 


1 
8 


4 


(2,2,0,0) 


(JVi - 1) ivf (iVi +1) (JV2 ~ 1) ]V| (]V2 + 1) 


{JV-1)JV^ 


1 




8{]V-2)(JV-1)]V(JV+1) 


(]V-2)(JV+2)(JV+3) 


8 




(2,1,1,0) 


(iVi-2){]Vi-l)JVi(]Vi+l)(JV2~2){]V2-l)JV2(iV2+l) 
8{]V-2)(JV-1)]V(JV+1) 


JV-1 
N+3 


1 
8 




(1,1,1,1) 


(iVi-3)(JVi-2){]Vi~l)]Vi(JV2~3){]V2-2)(JV2-l)]V2 
24{]V-3)(JV-2){]V-1)JV 


N-3 
N+1 


1 

24 




(5,0,0,0,0) 


JVi(iVi + l)(iVi+2){JVi+3){]Vi+4)JV2(]V2+l){]V2+2)(JV2+3)(JV2+4) 
120]V(JV+l)(JV+2)(JV+3)(JV+4) 


JV+4 
JV+9 


1 
120 




(4,1,0,0,0) 


(JVi -l)]Vi(JVi + l){JVi+2){]Vi+3)(JV2-l)Af2(-'V2 + l)(JV2+2){JV2+3) 


N+2 
N+7 


J_ 

30 




(3,2,0,0,0) 


(JVi~l)JvJ(]Vi+l){]Vi+2)(JV2-l)]v|{JV2 + l){iV2+2) 


(JV-l)JV(JV+2) 


J_ 








24 


5 


(3,1,1,0,0) 


(]Vi-2){]Vi-l)JVi(]Vi+l){]Vi+2)(JV2-2){]V2-l)JV2(Af2 + l)(iV2+2) 
20{]V-2)(JV-l)]V(JV+l)(JV+2) 


AT 
JV+5 


1 

20 




(2,2,1,0,0) 


(]Vi-2)(JVi-l)]v2(]Vi+l)(JV2-2){iV2-l)iV|(iV2 + l) 


(iV-2)(iV-l)JV 


1 




24(JV-2){]V-1)JV^ (iV+1) 


(JV-3)(JV+2)(JV+3) 


24 




(2,1,1,1,0) 


{iVi-3)(JVi-2){]Vi-l)JVi(]Vi+l)(JV2-3){iV2-2)(JV2-l)]V2(JV2 + l) 
30(JV-3){]V-2)(JV-1)]V(JV+1) 


JV-2 
JV+3 


1 

30 




(1,1,1,1,1) 


(jVi-4){]Vi-3)(JVi-2){]Vi-l)]Vi(JV2-4){]V2-3)(JV2-2){]V2-l)JV2 
120(A'-4){JV-3)(iV-2){iV-l)JV 


JV-4 
JV+1 


1 
120 



in (fT6l l to make this symmetry apparent 

„ n (A,+iVi-j)! iX,+N,-jy. 



X n 

i<j<j<;(A) 

iN--l{X)-iy 



N+l + Xi + \j-i~j 



n 



LA (A,+iV-Z(A)-z)! 



(17) 



We note that the obtained expressions for (sx) in terms of A^i 
and N2 yield zero if the length of the partition A is greater 
than n = mm{Ni,N2) so when averaging Schur function 
expansions one need not bother about the restriction? (A) < n. 

The Schur function average for f3 — 2 m terms of the chan- 
nel numbers A'^i and iV2, Eq. (fT4l i. has a simple structure, be- 
ing a ratio of polynomials 

(Sa)/3=2 = Ca 1 1 — (18) 



where the aj's are integers. Expression ( fTTI l for (sa)/3=i is 
less revealing. We found it useful to have the Schur function 
averages tabulated, see Table I for averages corresponding to 
partitions of to, m = 1, . . . , 5. This table suggests that 



(sa)/3=1 = Ca JJ^ 
J = l 



where the Uj are the same as in ( fTSl l and bj's are also integers. 
It would be generally desirable to understand the nature of the 
cancelations in ( fTTI l and to find a rule relating bj to A 4^ 



III. MOMENTS AND CUMULANTS OF THE 
CONDUCTANCE AND SHOT-NOISE 

We now apply the results obtained in the previous two sec- 
tions to calculate the moments of the conductance and shot- 
noise in a closed form. The final expressions involve summa- 
tion over all partitions of r in the case of conductance and 2?- 
in the case of shot-noise, with r being the order of the mo- 
ment. The cumulants Kr can be obtained from the moments 
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/ir with the help of the well-known recursion 

Kr = ^J.r _^.j^J■r~JK■j . (19) 

This method is well suited for analytic computations of lower 
order cumulants and also can be straightforwardly imple- 
mented in a computer algebra system for computations of 
higher order cumulants symbolically. Since a number of the 
partitions of the given r grows ~ exp{7ry^2r/3}/(4r\/3) 
only asymptotically at r oo (i.e. slower than pure exponen- 
tial), our method is very efficient for computing the cumulants 
up to reasonably large orders, as discussed below. 



The second sum on the right is over all partitions of r. From 
(ISTT l one easily obtains all moments of the conductance: 

(ffO =H ^ca(sa), r-1,2,...,. (22) 

Ahr 

For /3 = 2 this expression together with Eqs. (fT4l i and (fT5] l 
reproduces the recent result of Novaesi^l 

With Eq. ( |22] | in hand, one can obtain cumulants by making 
use of Eq. ( fT9] l. On this way we have successfully reproduced 
the first four cumulants which have been obtained before (ex- 
actly for any /3). For the reference purpose, we state explicitly 
the conductance variance'! 



A. Conductance 

The moments of the conductance can be obtained from the 
generating function 



(20) 



The desired Schur function expansion for the exponential 
function e^* can be read from ( fTTT i by choosing r, = 
there. Throughout this paper we use the convention that 
l/j- — for i < 0. The factorial determinant in (fTTT l can 
be evaluated by elementary transformations on its rows or 
columns and the answer turns out to be exactly the coefficient 
c\ introduced in ( fTSl l. Recalling that the Schur functions are 
homogeneous, thus sx{tT) — t^^^sx{T), one turns the Schur 
function expansion of the moment generating function into the 
following series in powers of t: 



(21) 



2{Ni + 1){N2 + 1) 
var(g) ^ I N{N +1){N + 3) ' 
(g) " ] 

{N- 1)N{N+1)' 



(5 = 2 



(23) 



with(5) =iViiV2/(iV+|- 
and the third cumulanti^ 



1) being the conductance average. 



4[1 - (jv^ _ N2f] 
iid')) J (A^-l)(A^ + l)(A^ + 5)' 
var(g) ] 2{Ni - iVa)^ 



(3 = 2 



(24) 



{N -2)NiN + 2)' 
which is a measure of the skewness of the probability distri- 
bution. An explicit expression for the fourth cumulant is quite 
lengthy for arbitrary iVi 2 and the corresponding large N ex- 
pansion can be found in Ref.'^ (see also below). However, 
in the particular case of Ni = N2 = n, it can be simplified 
further to the following compact form: 



Ahr 



3(4n4 + 20n3 + 43n2 + 53n + 24) 
((g^)) ^ J ~in+ l){2n - l){2n + l)2(2n + 3){2n + 5)(2n + 7) '' 
var(g) I 3 



2(2n - 3)(2n - l){2n + l)(2n + 3) ' 



/9=1 



(i = 2 



(25) 



Let us now discuss higher cumulants of the conductance, 
((g'")). Their explicit expressions are cumbersome and we 
consider mainly the physically interesting cases of the small 
or large channel numbers. In the quantum regime of a few 
open channels, we have found that these cumulants do not 
show a pronounced decay with increasing r, see Fig. 1 . In this 
case, the distribution function is strongly non-Gaussian. How- 
ever, as the number of channels in the both leads increases, the 
system approaches the semiclassical ('metalic') regime where 
one should expect"*'' the following dependence of the cumu- 



lants on the total number of channels, N = Ni + N2'- 

iif)) - {gf-- - N'-'. (26) 

The same scaling is generally applicable to any linear statis- 
tic on transmission eigenvalues (e.g., shot-noise), implying a 
Gaussian distribution in the limit — > cx;'*^-'''^ (see, however, 
the next section for discussion). 

We have performed the asymptotic analysis of our exact 
RMT expressions in the limit when both TVi 2 ^ 1. It sug- 
gests that the leading order term in the 1 /N expansion of the 
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FIG. 1 : The third to tenth cumulants of the conductance for chaotic 
cavities with a few open channels in the case of preserved (/3 — 1) or 
broken ((3 — 2) time-reversal symmetry. The corresponding values 
(NijNi) of the channel numbers are as indicated in the legend. 



that the Gaussian distribution is approached in this case much 
faster as compared to ( l26b - dZTl i. 

Generally, we note that in the symmetric case all odd cu- 
mulants at /? = 2 must vanish identically, as it follows by the 
simple symmetry argument'^ (indeed, the joint distribution (|2| 
becomes then symmetric under the change of all Tj 1^7} 
implying the symmetry of the conductance distribution about 
its mean §). It has been recently checked^'' that representation 



at /3 = 2 satisfies this property. For the even cumulants, 
the 1/n expansion of our exact expressions gives the follow- 
ing leading term at n 1 (fc > 2): 



(2fc-l)! 
4(4n) 



2k 



(28) 



that agrees with the recent result by Osipov and Kanzieper— 
obtained by a completely different method. In the case of /? = 
1, we have found with the help of symbolic computations in 
Mathematica that 



(r - 1)' 



1 , odd r 

-2n(r-3)!' 



(29) 



even r 



for r = 3, 4, . . . , 16. Correspondingly, we put forward the 
conjecture that Eq. ( |29] ) holds for all r > 3. 



r-th cumulant, r > 3, has the following general structure: 



(9) "(/3/2)' 



-G 



2(r-3): 



(27) 



where Gm{Ni, N2) is an independent of jS homogeneous 
symmetric polynomial of order m, see Table II for the first 
four ones. The expression on the rhs in dZTl i is of the order 
of N^^^', being in agreement with the above estimate ( l26l l ob- 
tained within a different approach (weak localisation diagram- 
matics). The next to leading order term in the 1/iV expansion 
of {{g^))/ (g), the so-called weak localization correction, is of 
the order of iV^''. It vanishes for systems with broken time- 
reversal symmetry {(3 ^ 2). Further terms in this expan- 
sion can be easily computed as well if necessary. 

In the special case of symmetric cavities, Ni — N2 — n, 
the leading term ( |27] ) in the 1 /N expansion of the cumulants 
vanishes for all r > 3 and so does the next-to-leading term 
of the expansion of any odd cummulant (independently of (3, 
it contains a factor (iVi — iV2)^ explicitly). This indicates 



TABLE II: The first four polynomial (iVi , 7V2), Eq. W\ . 



G2ir-3){Nl,N2) 



1 

- 4N1N2 + Ni 
Nf - lONfNi + 22N?Ni - 10iViiV| + Ni 
iVf - ISNfNi + 88NtNi - 150 N^Ni + {Ni ^ N2) 



B. Shot-noise 

Having an aim to find also the joint moments of the con- 
ductance and shot-noise, we consider the generating function 
for the moments of ag + p: 



Fit, a) 



{a+l)Ti+tTf 



(30) 



The moment generating function of shot-noise is then simply 
given by Fp {t) ~ F{t, 0) whereas that of the conductance fol- 
lows as Fg{t) — lima^oo F{t/a, a). At finite a, the quantity 
ag + p has a physical meaning of the total noise including 
both thermal and shot-noise contributions, with a being then 
the known function of the temperature and applied voltage.^ 

The exponential function in ( |30] | can be expanded in Schur 
functions s\{T) with the help of the general identity (fTTT l. On 
multiplying two exponential series, one obtains 



F(i,a) =^CA(i,a)(5A), 

A 

cx{t, a) = det{7rA,_i+j(t, 
where Trr{t, a) are polynomials in t and (a + 1), 



Lr/2J 

Trr{t,a) = ^ 

fc=0 



fc! {r-2k)l 



(31a) 
(31b) 

(32) 



In order to extract from this the moments of -I- p one 
needs to expand the coefficients c\ {t, a) in powers of t. After 
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some algebra, see Appendix lAl for details, one arrives at the 
desired expansion 

oo r 

i^(i,a)=^r^(-ir(l+a)'-™ 

r=0 m=0 Ahr+m 

(33) 

where 



E 



1 



}. (34) 



The determinant on the rhs ( [34] i can be evaluated in terms of 
the partition A leading to an explicit expression for the coef- 
ficients f\,rn, see Eq. (IA2l i. In the particular case of m = 0, 
/a,o is just the coefficient cx given by ( fTSl l. We note that 
depend only on A and m and not on n. The summation indices 
kj in (l34l l run over all integers from to m and are not subject 
to any ordering. From expansion (l33T l. one easily finds that 
the r-th moment of the total noise reads as follows: 



{(ag+pY) = rlJ2i-iy^l- 



-a 



m=0 
Ahr+m 



(35) 



where the second sum is over all partitions of r+m. The joint 
moment {g^p'^'^^) of the conductance and shot-noise is then 
given by Eq. (|35] |. with being replaced by the bino- 

mial coefficient ('^^.™). It is interesting to note that by setting 
a = — 1, one also obtains the moments of the sum of squares 
of the transmission coefficients: 



Ah2r 



To the best of our knowledge the above formulas have not 
been reported in the literature before. 

We now focus on the analysis of the shot-noise cumulants. 
Expansion ( l35l l successfully reproduces the general (3 results 
for the shot-noise average'' 



(3 var(5) 



(36) 



and for the shot-noise variance'^. The explicit expression for 
the later is rather lengthy (see Ref.'^ for the corresponding 
large N expansion) but turns out to be quite compact in the 
particular case of Ni — N2=n: 



U2n^ + 91^2 - 49n - 36) 



var(p) 



2(n+ l)(2n- l)(2n + l)(2n + 3)(2n + 5)(2n + 7) ' 

- 9n2 + 3 
An{2n - 3)(2n - l){2n + l){2n + 3) ' 

I 



/3 = 1, 
/3 = 2 



(37) 



Higher cumulants of shot-noise, ((p^)), similarly to those 
of conductance, are non-vanishing when the number of chan- 
nels is small, implying a strongly non-Gaussian distribution 
also in this case. In the opposite limit of the large number of 
channels, A'^i 2 ^ 1, we have found the leading term of the 
1 /N expansion to have the following structure: 



shot-noise (r > 3): 



(p) {l3/2y-^ 

x{N^ - 4:NiN2 + N^f 



(38) 



with Prn{Ni, N2) being an independent of /3 homogeneous 
polynomial of order m, see Table III. The next order term of 
the expansion has been found to have similar structure to that 
of the conductance (explicit expressions being, of course, dif- 
ferent), thus the same conclusions apply for this term, too. 

We consider now the particular case of Ni = N2 ^ n. 
In contrast to the conductance, both even and odd cumulants 
of shot-noise are non-vanishing at finite n, even for systems 
with broken time reversal symmetry {fi ~ 2). In the limit of 
n ^ 1, the 1/n expansion of our exact expressions suggests 
the following asymptotic behavior of the r-th cumulant of the 



0=2 



and 



(r-l)l J 1' 

r-II ' 



(39) 



(40) 



even r 



Til 



TABLE III: The first four polynomial Pm{Ni, N2), Eq. lUll. 

r P4(r-3){Nl,N2) 

3 1 

4 iVi* - IGN'iNi + 34NfNl - WNiNi + iV| 

5 Aff - SSNlNi + 385Af A| - 1344Af A| + 2008Ai*A2* 

+ {Ni ^ N2) 

6 Nl^ - 66A1" A2 + 1345Ai"'A| - 11680Af A| + 49699Af A| 

-114598Ai^A| + 150662Af A| + (Ai ^ A2) 
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FIG. 2: The conductance distribution | |4U for chaotic cavities with 
preserved (/3 = 1) or broken (/3 = 2) time-reversal symmetry. The 
number of channels is fixed to A'^i =2 in one lead and varied in 
the other, = 2, 3, 4 or 5 (solid, dashed, dotted or dash-dotted 
lines, respectively). Non-analyticity of the distribution (a cusp point 
at g = 1) becomes less pronounced as increases. 



Equation ( l39b agrees with the very recent result^* obtained by 
a different method. We have been able to verify by symbolic 
computations in Mathematica that Eq. ( i40l i holds up to the 8- 
th cumulant and, thus, conjecture it to hold for any r > 3. 



IV. DISTRIBUTION FUNCTIONS 



ductance, Pg^\x) 



We consider now the distribution function of the con- 
— (S{x — g)), and that of shot-noise, 
p y^) — {S{x — p)), with g and p being defined in ([T]). 
Explicit expressions for the conductance distribution can be 
found in the particular cases of n = 1 , 2. At TVi = 1 and N2 = 
K > 1, Eq. ^ readily gives P^'^^^^ig) = {(3K /2)gP^/^-'^ 
for < g < 1, and zero otherwise. In the case of iVi = 2 and 
arbitrary N2 — K > 2, the conductance distribution can also 
be found by performing integrations that feature in the defini- 
tion (see (IB4b in AppendixlBl). with the final result being 

PS=2i9) - Kgf^^-'[X,-i-l)(f^^-^y^X2e{g-l) 



^i2(^-)^^-9iliK-^)+J^^-f^K)] (41) 

for < g < 2 and zero otherwise, see Fig. 2. Here 
Bz{a, b) is the incomplete beta-function, 8(x) stands for the 
Heaviside step function, and the constants Xi 2 are given by 

y _ r[f}{K+l)/2+l]r(pK/2) , y _ r[/3(jf+l)/2+l] p 

^1 ~ r(/3/2)r(/3/f) ^2 — r(/3)r[0(A:-i)/2] ■ 

pression dJU holds for arbitrary positive integer /3. In the par- 
ticular case of (3 — 1, it can be simplified further, yielding 
P^,IL2{9) = ^K{K+1)[{^)^-^ - (.g-l)(^-i)/2e(.g-l)], 
in agreement with RefM. 

It is possible to find explicit expressions for the distribu- 
tion function of the conductance beyond the cases discussed 
above. However, the final answers become more cumbersome 
as the channel numbers grow, thus being almost of little prac- 
tical use. For the shot-noise distribution, the situation is not 
satisfactory even for small channel numbers: we are not aware 
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FIG. 3: Comparison of the exact conductance distribution (Eq. | |4U . 
black solid lines) with the Edgeworth approximation (Eq. l |43l l, red 
dashed lines) for chaotic cavities with preserved (/3 — 1) or broken 
(/3 = 2) time-reversal symmetry. The number of channels is fixed to 
A^i = 2 in one lead and varied in the other, N2 — 3, 4, 5, 7 (lines 
from left to right, respectively). The purely Gaussian approximation 
would produce much stronger systematic deviations, thus being not 
sufficient at small channel numbers at all. 



of explicit results for the shot-noise distribution except for the 
simplest case of A^i ^ — 1.^^. One practical way to solve 
this problem is to construct approximations to the distribution 
functions in terms of cumulants by making use of the Edge- 
worth expansion. It turns out that such approximations are 
fairly accurate in the bulk of distribution even for small chan- 
nel numbers. 

In the limit A^i 2 ^ 1 the conductance and shot-noise dis- 
tributions follow the Gaussian law 



(t)o{x) 



1 



(42) 



where /i and are the corresponding mean value and vari- 
ance, respectively. The Edgeworth expansion is a ex- 
pansion around the Gaussian law. Denoting the r-th cumu- 
lant by the first correction to the Gaussian law is given 



by 01(2;) 



-^Ksd'^cpoix) and the next one is given by 



'f>2ix) = (jfK4(9^ + ^K^d^)(f>o{x). Higher order connections 
involve higher order cumulants^ ^ Restricting ourselves to the 
first four cumulants, we get the following approximation to 
the distribution functions of interest: 



(43) 



k=0 



The advantage of the Edgworth expansion is that it is a true 
asymptotic series, with the controlled eiTor (e.g., in our case 
the error of approximation (|43T l is estimated to be of the order 
of ■^). With the higher order cumulants being readily avail- 
able from Eqs. (i22] i. ( |35] | and ( [T9l ), one can easily improve the 
accuracy of the approximation by adding higher order correc- 
tions if needed. 

We found it instructive to compare the above approximation 
(l43T l with the exact results which are available for the conduc- 
tance distribution. To our surprise, we found the Edgeworth 
approximation to be fairly accurate already for Ni = 2 and 
N2 > 4, see Fig. [3] The agreement between the approximate 
and exact distributions gets even better for iVi.2 > 3. 
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The Edgeworth approximation fails near the edges and also 
if the distribution has strong singularities. Therefore, it is 
tempting to look for an alternative and exact representation for 
the distribution functions of the conductance and shot-noise. 
We note that at any finite number of channels each of these 
distributions has a finite support, namely, < g < n = Ig and 
< p < J = Ip, being identically zero outside this region. It 
is, therefore, natural to represent the distribution functions in 
the following form (henceforth, variable and index x = g^p): 



(44) 



as the Fourier series over the interval of support, cf. Schur 
function expansion (IT). 

We now show that the Fourier coefficients are given by 
Pfaffians'*'' as follows: 

(m) = Im Pfaff [^i'^) (m)] . (45) 
■Np 

Before establishing explicit forms of the anti-symmetric ma- 
trices a'£'^ (m), it is useful to note the following. The Fourier 

coefficients turn out to decay generally as cjf'' (m) ~ m""^^ 
as m oo, where exponent > 1 depends on the case con- 
sidered. This readily leads to the observation that the every- 
where continuous distribution function (l44l i contains (rather 
weak) singularities at integer points of division of the support 
interval [0, Ix], as certain derivatives become discontinuous at 
these points (see further Ref. ' ' for an alternative geometric in- 
terpretation of these singularities). Such a non-analyticity is 
expected to become less and less pronounced when the num- 
ber of channels grows, as the bulk of the distribution is de- 
scribed then by a Gaussian law. However, the asymptotic be- 
havior near the edges of support is always characterized by a 
power lawii. Consequently, even at n 1 the distribution 
remains to be weakly singular at the junction of Gaussian and 
power-law regimes^^, see also Ref.^^. 

We start the derivation of (l45T l with the simplest case of 
unitary symmetry, (3 = 2, and consider first the conductance 
distribution 

P^'Hg) = j^Jd[T]l[Tr'A{Trs[g~Y.T. 

j=i 1=1 

Writing A(r) = det{Tr^}, with i, j = 1, . . . , n, and sub- 
stituting the Fourier representation of the (5-function, 6{g — 
T.-) = / |^e-*"9+*" 5:. Qjjg ^^jj interchange the order 
of integrations and then apply the Andrejeff identity to per- 
form the integrations over T/s. This readily yields 

I /"OO 7 

P^'Kg) = ^J ^e-^-^det[AfK^)], (46a) 



kl 



^rprj-ia-^k-\-l — 3^iujT 



(46b) 



Since we know that the conductance distribution has a support 

(2) 

only in [0, n] it is actually more convenient to expand Pg ' (g) 



in functions ^j2/ns\n{raT:g/n), with to = 1,2, . . ., which 
form a complete and orthonormalized set on this interval. Fur- 
thermore, due to factorization of the exponential function, one 
has to express first everything in terms of exp{irmrTj / n) and 
then take the imaginary part. As a result, we arrive at the final 
answer cast in the form of (05]), where 



Pfaff[A(2)(-^)] 



det[i(2)(TO)] 



(47) 



and if ^ (to) EE A^^\^) is the discrete analogue of ( l46b] i. 

In the case of orthogonal symmetry, (3 = 1, the derivation 
goes along the same lines but has to be done separately for 
even or odd n. In the case of even n, we first consider the 
Tj's in a special order Ti < T2 < . . . < Tn (hence nl). 
Then it is useful to represent the Vandermonde determinant 
as a Gaussian integral over two kinds of Grassmann variables 
and integrate out further one set of them using the method 
of alternating variables^, see also Ref. for relevant details. 
The resulting expression acquires then the symplectic (anti- 
symmetric) structure automatically and the remaining average 
yields and (gill, with 



^«(to) 



dTi / dT2sign(r2-Ti)T^ 

JQ 

^rpa+l~2^i(m7T/n){Ti+T2) 



a+k-2 



(48) 



In the case of odd n, one has to increase artificially the number 
of the Grasmannians by one and proceed as before, with the 
final result being 



AW(to) 



odd n 
kl 



Bk,n+1 



(49) 



Here the nxn matrix A'^^^ is given by Eq. ( 1481 ), and the n- 
dimensional vector B is 

Bk,n+1 = fdTT'^+k-^^^im.MT ^ ^^^^^^ 

Jo 

For the sake of completeness, we also state the result in the 
case of symplectic symmetry, /3 = 4: 



kl 



(l-k) dTT 



na+k+l- 



■4gi(m7r/n)T 



(51) 



where now k,l = 1, . . . , 2n. It can be obtained by reduc- 
ing the fourth power of the Vandermonde determinant to the 
calculation of a Vandermonde determinant with the doubled 
dimension 2n subject to the additional (5-function constraints 
for the corresponding pairs of eigenvalues. The rest is as in 
the orthogonal case above. 

We note that it turns out to be possible to find explicit ex- 
pressions for the conductance distribution in the case ofn = 3 
and 4 by evaluating the above Pfaffians analytically and per- 
forming the corresponding Fourier transformation exactly. 

The expressions derived in this section are also well suited 
for numerical calculations. It is important to note in this re- 
spect that, as the Pfaffian is defined as an analytic square 
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root of the determinant of an antysimmetric matrix A, a spe- 
cial care has to be taken to decide for the complex determi- 
nant which sign has to be chosen. To overcome this dif- 
ficulty, we outline the following general procedure. First, 
we multiply the matrix A with the symplectic unit Z = 
diag{— 1(72, • ■ • , — io'2}, (72 being the Pauli matrix. Then we 
note that the matrix ZA is selfdual, ZA — AZ, implying 
that its eigenvalues come in pairs. Taking the product of all 
eigenvalues of ZA only once (this is exactly the quaternion 
determinant of ZA), we obtain finally the Pfaffian of A. As 
an illustration of this procedure, let us consider the simplest 
example of calculating the Pfaffian of o ) ^ i(72a which 
is equal to a. The matrix (— icr2)(i(T2a) = diagja, a} has 
obviously eigenvalues (a, a), so that taking a once yields the 
Pfaffian a. This is exactly the way how Pfaffians can be easily 
computed numerically. 

Generally, one can derive the distribution function of any 
linear statistic onT, x — J^j fx{Tj) with a given fx{T), in 
complete analogy with the above lines. It yields representa- 
tion ( |44| ). where matrices Ai^'^ (to) are given by the above ex- 
pressions (|47]| - (BTI) in which all the factors exp{i{rmT/n)T} 
have to be obviously substituted with cxp{i{'imT / Ix) fx{T)} , 
with Ix being the length of the corresponding support interval. 
In particular, for the case of shot-noise it amounts to substitut- 
ing there with exp{i(4TO7r/n)r(l— T)}. 



V. ASYMPTOTICS 

The asymptotic behavior of the distribution functions near 
the edges is characterized by a power law dependence 



(52) 



where both the exponents and the pre-factors can be deter- 
mined exactly at arbitrary iVi 2 and any /3, as shown below. 

In the case of the conductance, the exponents £g and Vg have 
been already reported previously ' ' , being given by 



ig = an + ^{n — l)n — 1 
rg = (n-l)(l + fn) 



(53) 




FIG. 4: The distribution of the conductance {Pg{g), left) and that of 
shot-noise (Pp{p), right) for chaotic cavities with preserved (/3 — 1) 
or broken (f3 = 2) time-reversal symmetry. The number of channels 
are Ni = 2 and ~ 4. Dashed lines correspond to the exact 
distributions whereas solid lines show the corresponding asymptotic 
behavior near the edges, see the text for details. 



The behavior of the distribution near the right edge can be 
analyzed in a similar way, yielding ( |52] | with 



Rg - 



1 



r 



Ji-1 

n 



r[l + a + f(n+j-l 



(55) 



In the case of shot-noise, the corresponding exponents are 
found to be as follows: 



£p = an + ^[{n - If - a] - 1 
rp = | + |(n-l)n-l 



(56) 



where cr = 1 or for n even or odd, respectively. The above 
expression for Vg was already known" whereas that for £p is 
new, being in agreement with the recent large n result. The 
corresponding constants can also be found exactly. We refer 
to Appendix IbI for further details, including the discussion of 
correction factors to the leading asymptotics given above. 

As an illustration of the obtained results. Figure 4 shows a 
comparison between the exact and asymptotic behavior of the 
conductance and shot-noise distributions in the case of chaotic 
cavities with Ni = 2 and N2 = 4: channels. 



To determine the constant Lg, we consider Pg{g) at 5 < 1. 
In this case, the upper limit of the integrations over T^'s in 
^9 (.9) ~ {^{9 ^ J2i T^i)) i^^y be replaced with g (due to the 
(5-function). Scaling further all Ti gTi and calculating the 
powers of g there, one readily gets Pg{g) = Lgg^^^, with Lg = 
((5(1 — ^ Ti)). It is worth emphasizing that this is the exact 
expression of the conductance distribution at < g < 1.— 
The integral for Lg can be then calculated by the standard 
RMT methods, see Appendix iBl and is given by 



L„ 



1 



n 

J=0 



r[i 



f(n+j-l) 



r(i 



(54) 



VI. CONCLUSIONS 

In this work, we have presented a systematic study of statis- 
tics of the conductance and shot-noise in chaotic cavities. Our 
approach is based on expanding symmetric functions in the 
transmission eigenvalues, of which the conductance and shot- 
noise are two examples, in Schur functions and then applying 
a generalization of Selberg's integral due to Hua to evaluate 
the averages. This leads to explicit formulas for the conduc- 
tance and shot-noise cumulants in terms of the channel num- 
bers for chaotic cavities with preserved (/? = 1) or broken 
{(3 = 2) time-reversal symmetry. 
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For lower order cumulants our formulas reproduce the pre- 
viously known exact results. We have performed an asymp- 
totic analysis of the cumulants in the regime when the number 
of channels in both leads increases. It suggests that, generi- 
cally, the r-th cumulant decays as l/N"^^"^ for r > 3 with the 
leading term containing (A^i — A^2)^ as a factor; see Eqs. (|27] | 
and (|38] l. This implies that the convergence to the limit- 
ing Gaussian law is faster in the case of symmetric cavities, 
Ni — N2 — n. In this case, we have been able to analyze 
the first 16 cumulants of the conductance and the first 8 cu- 
mulants of shot-noise in the limit n ^ 1, obtaining the lead- 
ing order term in the 1/n expansion explicitly in terms of the 
cumulant order, r, and the channel number, n. For the sys- 
tems with broken time-reversal symmetry, the r-th cumulant 
decays as for both conductance and shot-noise, with all 
the odd cumulants of the conductance being identically zero. 
Our results in this case, Eqs. (l28T l and (|39] l, agree with those of 
Osipov and Kanzieper^^'^^ obtained recently by a completely 
different method. For systems with preserved time-reversal 
symmetry, we have found that the r-th cumulant decays as 
1/n'' for odd r and 1/n''^^ for even r, r > 3, for both the 
conductance and shot-noise. This staircase effect in the rate 
of the cumulant decay seems to be a novel feature which has 
not been reported in the literature before. One of its apparent 
consequences is that the convergence to the limiting Gaussian 
law is slower for systems with preserved time-reversal sym- 
metry. We have also put forward our explicit formulas ( |29] l 
and (l40l l for the higher order cumulants in the whole range of 
r > 3 as a conjecture. Proving this conjecture seems to us an 



interesting and challenging open problem. 

As mentioned above, in the limit when the number of open 
channels in both leads increases, the conductance and shot- 
noise distributions are described by the Gaussian law. With 
higher order cumulants in hand, one can easily obtain next 
order corrections to the Gaussian law by making use of the 
Edgeworth expansion. We have found that such approxima- 
tions to the distribution function are fairly accurate in the bulk 
even for small channel numbers. We have also obtained an 
alternative and exact representation for the distribution func- 
tions, in terms of Pfaffians, which is suitable in the whole 
range of support, including the edges where the distributions 
have a power-law dependence. Such an asymptotic behavior 
have been investigated in detail, the powers and corresponding 
pre-factors being determined exactly at any (3 and A^i 2- 
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APPENDIX A: SCHUR FUNCTION EXPANSION, EQ. Ofi) 



In this appendix we expand the coefficients c\(t, a) — det {Tr\. ^i^j{t, a)} in powers of t. From ( |32] ) and the definition of 
determinant, 

" ( T^kifT , \Xi-i+a(i)-2ki j.\i-i+a(i)-ki 

det{..._.+,(t,a)}=;^(-l)HnE-^ ^ ^ 



n 



jr 



where K = J^i and the first sum is over all permutations a in the symmetric group Sn- Changing the order of summations 
again, one can fold the sum over permutations into a determinant again. This yields 



L|A|/2J , 

.(M)= E (-ir(l + E det _ . (Al) 



LIAI/2J 

m=0 fei+...+fc„=m K'^i-V^i " ' J ^"-J^Jij^l 

The second sum here is exactly the coefficient /A,m that appears in Eq. (|30] |. The determinant in jAlb can be evaluated by 
elementary transformations on its rows and columns, 

^^^j 1 \" _ ni<,<,<„(A.-*-A,+j-2fc, + 2fc,) 



fc,!(A, + 2k,)\ J n;=i kj\ {n + A, - j - 2kj)\ 

Note that for any n > 1{X), where /(A) is the length of A, one can safely replace n in the above expressions by /(A). Hence, the 
sum over the n-tuples (fci, . . . , fc„) of integers in (lAll i can be replaced by the sum over the Z(A)-tuples (fci, . . . , fc;(A)), yielding 

/A,m = > TTT^ ■ (A/) 

n -LI k,l (/(A) + A, - J - 2k,)\ 
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We now substitute the obtained expression for c\ {t, a) back in dSTI ) to obtain 

oc. L'72J 

Fit, a) = J2J2 + ay-'-^f-"^ hMp, 

r— m— Ahr 

hence the desired Schur function expansion dSOb follows after changing the order of summations. 

APPENDIX B: ASYMPTOTIC CONSTANTS AND CORRECTION FACTORS 
Conductance, left edge: To calculate Lg = ((5(1 — Tj)), we use the following result due to Mehta^ (see p. 361): 

f TTixx-^lAfx)!^ e(i-i:.a:0 r(7) " r[. + |(n-.)]r(i + fo 

o-Lr^ ^' ' ^ %i-E.xO^-~r[, + an+f.(n-l)]fJ r(l + f ) 



(Bl) 



Making then use of the known identity lim^^o+ 7(1 ~ ^©(1 — x)Q{x) — 5{\ ~ x) and noting that 7r(7) = r(l + 7) ^ 1 
at 7 ^ 0, one can readily get from ( IB lb the following result (valid at any (3 > 0): 

1 'fT^r(a + fj)r[i + |(j+2)]^ 



■■■ d[T] n A(T)i^,5(i -Y^T^)^ — - n 



r[a7i+fn(n-l)]/^^Q r(l + f ) 

After dividing (IB2b with the normalization constant (|3]l, this yields Eq. ( |54| | of the main text. 

Conductance, right edge: For considering the limit g ^ n from below, we make a transformation Ti — > 1— T!; in the integral 
for Pg*''^ (g) M-^j d[T] Hj (1 - T^)""^ | A(r) |''(5(n - g - T^) and then consider small. Now one sees that in the region 
n > 5 > n — 1 the upper bound of the integral can be replaced by n — 5. Further rescahng Ti (n — g)Ti yields 

Pj'^)(5) ~i?g(n-g)("-i)(i+^")e-("~i'("-s), (B3) 

where Rg is then given by Eq. (|55] ) (use result ( IB2I ) at a = 1 and divide it by Afp). The correction factor e^'^"^^'^"^^) comes 
from expanding Ylii^ ~ {n^9)Ti)°'^^ linearly in the exponent. Thus we have systematically expanded the log of the positive 
quantity Pg^\g) for small [n — g) including the term of order [n — g). Corrections are of relative order [n — g)^. 

Conductance, 1 < g < 2: Since we know Pg^\g) in the interval < 5 < 1 exactly, one can consider asymptotics in the 
interval 1 < g < 2 for g — > 1. To this end, we first arrange the integration variables as Ti < T2 < • • • < T„, then rescale 
Tfe TnTk for k = 1,2,..., n— 1, and finally perform the T„ integration. As a result, one arrives at the following exact 
representation of the conductance distribution at < g < 2: 

PfKg) = _ P{g)l (B4) 

where F(.g) = {n/Np) J d[T]\A{T)f Ht^ Tr'a - T^^i^ + E^i' T,) -'^-(-^W^Q^g 1 _ ^-1 j^^) . :„ the par- 
ticular case of n = 2, this integral is one-dimensional and can be easily evaluated explicitly, resulting in (|4TI) . At any n > 2, 
P{g) can be handled with the same methods as before, yielding the following asymptotic behavior at small g — 1: 

P{g) ~ e(g - l)H{g - i)"("-i)+l(«-2)(«-i)e-^(9-i) ^ (B5) 

where the constants H and J are given by 

n[a + §{n-l)] T[l + a + Un + j - I)] 



H 



r[l + a(n - 1) + |(n - 2){n - 1)] f^^ T[l + f (j + 1)] 



rn + ^(7 + i)i 



and J = (" "'"''^"_|!^"„^)^_^^^'L^)|,f^^^^^ . These are the expressions used to make plots on Fig. 4. 

Shot-noise, right edge: Shifting all T, ^ i + T, gives P^^\p) = Af^^J d[T] Uji^+Tjr'^\MT)\f>S{'i - JJ - T^), 
where the integration is now over the n-dimensional cube centered at origin: — ^ < < |. Considering the right edge of the 
support of Pp'^-' (p), thus J — p ^ 0, one can then rescale Ti Ti — p and obtain: 

^ Pp(n/4 -p)"/'+"'"-')''/'-'e-("/4-P)^ (B7) 
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where Rp = (2(i-")"/AA^) J^^- ■ -f^^ d[T]\A{T)f5{l - J^i Tf). This integral can be further reduced to a Gaussian type 
integral found in Mehta^ (see p. 354). After some algebra, we finally arrive at 



2("-i)"+/3(«-iW4r(rp + 1) 11 r(i + fj)r(a + |j) ' ^ ^ 



where is given by (|56] |. The correction factor e ^"^^^ with F = 2(a — 1) — 2(a — + §(n — 1)], is found as before 

by expanding the rest in powers of ^ — p and retaining the linear part in the exponent. 

Shot-noise, left edge: Forp 0, one gets positive contributions to Pp^\p) from all 2" corners of the integration cube. As 
function T{l — T) is not monotonous, it is more convenient to treat the contributions from all the corners separately. To this end, 
we make for {n — m) variables T^. the transformation ^ (1 — T/.) and then let all Ti run from to i. This gives 

P m=0 ^"^^ 1=1 k=m+l 

X n \T^-Tjf n \n-Tif6{p-J2ni-T.)). 

l<z<j<m m+l</c</<n i 

The factor (^) appears since one has (^) equivalent corners, and 2^ = X!m^o (m)- ScaHng Ti ~> pTj yields 

P(«(p) QQ(m)p^(™)-ie^^(™), (B9) 

m=0 ^ ^ 

with 7(to) — n + m(a — 1) + ^[m{m — 1) + (n — m)(n — m — 1)] and 

r[7(m)]li r(a + jf) ii r[l + |(n-m + j)]r[l + f(n-m + j + l)]' 

2m /3 
^;(to) = 3 - a + f3(n~2m~ 1) + ——[a- 1 + l3(2m- n)]\a+ -(m- 1)]. (Bll) 

7(m) 2 

For the given even n = 2fc or odd n = 2k + 1, k > 1, the function 7(to) has a minimum at m = fc that gives the leading 
exponent i!p = -f{k) — 1 stated in Eq. ( |56] |, and = (^)Q{k). 
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